Majorana edge states in two atomic wires coupled by pair-hopping 



Christina V. Kraus,^'^ Marcello Dalmonte,^ Mikhail A. Baranov/'^'^ Andreas M. LauchU,^ and P. Zoher^'^ 

^ Institute for Quantum Optics and Quantum Information of the Austrian Academy of Sciences, A-6020 Innsbruck, Austria 
^Institute for Theoretical Physics, Innsbruck University, A-6020 Innsbruck, Austria 
^RRC "Kurchatov Institute", Kurchatov Square 1, 123182, Moscow, Russia 
(Dated: February 5, 2013) 

We present evidence for the existence of Majorana edge states in a number conserving theory 
describing a system of spinless fermions on two wires that are coupled by a pair hopping. Our 
analysis is based on the combination of a qualitative low energy approach and numerical techniques 
using the Density Matrix Renormalization Group. We also discuss an experimental realization of 
pair-hopping interactions in cold atom gases confined in optical lattices, and its possible alternative 
applications to quantum simulation. 
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At present there is significant interest in identify- 
ing physical setups where Majorana fermions (MFs) [1 
emerge as a collective phenomenon in many-body quan- 
tum systems [2 . The motivation behind this search is 
two-fold: First, the existence of MFs is intimately linked 
to the concept of topological phases and their explo- 
ration. Second, MFs provide due to their topological na- 
ture a promising platform for topological quantum com- 
puting and quantum memory [3-5 . In a seminal paper 
Kitaev pointed out a route towards the realization of MFs 
in a simple many-body system [6]: A ID wire of spin- 
less fermions with a p-wave pairing can exhibit a topo- 
logically ordered phase with zero-energy Majorana edge 
modes. The key ingredient here is the coupling of the 
wire to a superconducting reservoir in a grand canonical 
setting, which is induced in complex solid state struc- 
tures via the so called proximity effect. Building on this 
result, a remarkable theoretical and experimental effort 
has been devoted in search of alternative settings sup- 
porting topological superconductivity in ID condensed 
matter systems, such as the combination of spin-orbit 
coupling, magnetic fields and s-wave interactions [7HT6]. 
Alternatively, Majorana physics can be observed with 
ID quantum gases coupled to a particle reservoir rep- 
resented by molecular condensates, taking advantage of 
the unique tools for control and measurements in atomic 
systems pTHTO] . 

In contrast, we propose and investigate in the present 
Letter an alternative approach to create Majorana edge 
states in a purely number-conserving setting [2QH22] . We 
consider the conceptually remarkably simple system of 
spinless fermions in two wires with single-particle in- 
trawire hopping, which are coupled via an interwire pair 
hopping (c.f. Fig. [l^). On an intuitive level the relation 
to Kitaev's model and existence of Majorana edge states 
is apparent, when we consider one of the wires as an effec- 
tive particle reservoir for the second wire. The essential 
element in our system is pair hopping between the wires, 
which breaks the U{1) symmetry associated with the con- 
servation of the particle-number difference between the 
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FIG. 1. a) Ladder Hamiltonian: Atoms in the a- and b- wires 
can tunnel individually along the x-directions, and can hop 
in pairs between the wires, b-e) Implementation of the pair 
hopping: b) The single wire is realized as a bipartite lattice 
of t and ^ fermions with Raman assisted tunneling (Rabi 
frequency Q). c) Ladder scheme as a combination of two 
wires with opposite energy off-sets. The dashed box denotes 
a single plaquette, with site indices indicated in parenthesis 
tam,tbm are the tunneling amplitudes from the a- resp. 
and 6-wire to the central sites. Atoms in the center (C) of 
the plaquette interact with strength U (shaded areas), d-e) 
Energy off-sets along the diagonal of the plaquette in for 
the l- resp. t species, and corresponding virtual processes 
indicating pair tunneling (see text), f) Time evolution of 
the state a\ ^^a^ ^^0) according to the microscopic dynamics 
in units t/tam (see text): the blue (red) curve indicate the 
pair population Pa,b{t) in the upper (a)/ lower (b) wire as a 
function of time. 



two wires, down to the Z2 parity symmetry, an ingredi- 
ent known to be crucial for the emergence of Majorana 
modes in the grand-canonical scenario. The purpose of 
this work is two-fold. First, we provide evidence for Ma- 
jorana edge states and related topological order using 
both field-theoretical arguments and detailed a Density- 
Matrix-Renormalization Group (DMRG) study j23l[2i]. 
Second, we show that the present setup with pair hop- 
ping has a natural implementation with cold atoms in 
state-dependent optical lattices [25 combined with Ra- 



2 



man assisted tunneling processes. 

The emergence of Majorana edge states in the super- 
fluid phase of the system is demonstrated on the ba- 
sis of the following criteria: (i) two degenerate ground 
states with different parities for the individual wires in 
the case of open boundary conditions (OBC), (ii) non- 
local fermionic correlations between the edges, coming 
along with (iii) topological order indicated by a degen- 
erate entanglement spectrum, and (iv) robustness of the 
above properties against static disorder. We also show 
that properties (i)-(iii) survive in the presence of a weak 
single-particle hopping between the wires, also support- 
ing the topological origin of the state. For experimental 
realizations, in particular with atoms, the last property 
could be crucial, as it shows the robustness against the 
most probable major imperfection. 

Model. We consider the following Hamiltonian: 



(1) 



where aj(a^), bj{b^j) are fermionic annihilation (creation) 
operators defined on two distinct wires a and 6, respec- 
tively, the first line describes intrawire single-particle 
hopping with the corresponding amplitudes ta^b (in the 
following we consider = ^6 = t as a weak asymmetry of 
ta^b does not affect the results qualitatively), and the last 
term is the interwire pair hopping with the amplitude W. 
The choice of the Hamiltonian (1), motivated by previ- 
ous considerations of the number-conserving setting [20] , 
stems from global symmetries and corresponding con- 
served quantities: Beside the total number of particles, 
N = Na^Nb = a]aj^b]bj, associated with the U{1) 
symmetry, there is another conserved charge - the parity 
Pi of one of the wires [say, the wire a. Pi = Pa = (~1)^''] 
associated with a Z2 symmetry [44] . The conservation is 
guaranteed by the last term in H allowing only hopping 
of particles between the wires in pairs, and is the key 
requirement to access a topological phase with MFs. 

Before presenting the analytical and numerical analy- 
sis of the Hamiltonian (1), let us give an intuitive picture 
based on the simplest system supporting fermionic Majo- 
rana edge states - the ID Kitaev quantum wire [6^ with p- 
wave pairing described by a mean- field BCS-like Hamilto- 
nian resulting from the coupling to a reservoir of Cooper 
pairs (see Ref. [6j for details). In our case, one could 
view one wire as a reservoir of pairs for the other wire 
and vice versa, and decompose the pair-hopping term 
in a mean-field manner as W Y-{a\a\j^-J)ibij^i + h.c.) 
Y,.{/\ba\a\^^- /^Ibibi+i^h.c.), where = W {aiai+i) 
and A5 = W {bibij^i) are non-zero pairing amplitudes 
which can be found by applying the standard Bogolyubov 
procedure. With this decomposition, the Hamiltonian 
(1) describes two Kitaev wires [6 , each of them having 
doubly-degenerate ground states with different fermionic 
parities Pa^b = =t for the a- and 6-wire, respectively, and 



carrying two Majorana operators corresponding to the 
edge-modes. Therefore, the ground state (GS) of the 
double- wire system (1) with a fixed parity Pi = Pa and 
a total parity P = (—1)^ is doubly degenerate. The two 
ground states can be connected by the product of two 
Majorana operators - one from each wire. Strictly speak- 
ing, long- wavelength fluctuations destroy long-range or- 
der in ID breaking the mean-field description even at 
zero temperature. In the considered case, however, this 
does not change the picture qualitatively (see Ref.[21J). 

Low- energy theory. Effective field theories based on 
bosonization [29l |30] represent a remarkable tool to in- 
vestigate the emergence of topological states and MFs 
in strongly correlated systems jll] [26H28] , and has been 
applied recently to number conserving settings [20l [21] . 
Here, we employ this formalism to qualitatively ana- 
lyze the low-energy properties of the Hamiltonian ([T]). 
We start with applying standard bosonization formu- 
las to introduce effective low-energy phase and density 
fluctuation fields (f^^'d^^ respectively, for each species 
J = a^b [45 . After introducing symmetric and anti- 
symmetric combinations, (fs/A = {^a i ^b)/V^j and ne- 
glecting contributions with high scaling dimensions, the 
bosonized Hamiltonian decouples into symmetric and an- 
tisymmetric sectors. The symmetric sector describes col- 
lective density-wave excitations, and is well-captured by 
a Tomonaga-Luttinger liquid Hamiltonian: 
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whilst the antisymmetric one is described by a sine- 
Gordon Hamiltonian I29|: 
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where Ko, and are the Luttinger parameter and the 
sound velocity, respectively, for each sector a = (A^S)^ 
and w (X W results from the pair hopping. It can be 
shown that the parity symmetry Z2 and the number con- 
servation are exactly retained at low energies in the an- 
tisymmetric and symmetric sector, respectively [201 [29]. 
We now discuss the qualitative phase diagram of the sys- 
tem by using standard Renormalization Group (RG) scal- 
ing arguments [29H3T] . Away from the strong coupling 
limit W ^ t (where terms with higher scaling dimensions 
may become relevant), the two sectors remain decoupled, 
so that one can analyze them separately. While the sym- 
metric sector is simply a theory of free bosons, the anti- 
symmetric sector displays richer physics, as it undergoes 
a phase transition from a gapless phase at = to a 
gapped, superconducting phase for 1^ > 0. In analogy 
with the continuum model of Ref. [20 , Eq. (|3| can be 
exactly mapped to the continuum version of the Kitaev 
wire [20 at the Luther-Emery point Ka = 2. As a result, 
the system with OBC displays a two- fold ground state de- 
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generacy, where the two states have opposite parities Pi , 
and support MFs at the boundaries [45]. Moreover, the 
single-particle correlation functions show exponential de- 
cay (alai-^x) — e"*^!^! in the bulk, signaling the presence 
of a finite superconducting gap. Away from the Luther- 
Emery point, the MF wave function overlap increases de- 
pending on {Ka — 2), the corresponding splitting in the 
GS degeneracy at finite system size being e~'^^, even- 
tually turning to power-law in the presence of certain 
kinds of perturbations [201422) . On the other hand, in 
the strong coupling limit \W\ ^ the presence of addi- 
tional terms (with higher scaling dimension) of the form 
w cos[\^47ri9A]{dx(ps)'^ — —'^{dx^s)'^ leads to a reduction 
of the sound velocity vs, resulting in phase separation. 

Numerical results Employing this low-energy picture 
as a guide, we now present a quantitative numerical in- 
vestigation of the Hamiltonian Eq. ([T]). We start with a 
brief description of the phase diagram of the system, and 
then discuss the criteria (i)-(iv) relevant for the existence 
of MFs. In the following, we set t = 1 as the energy scale. 

The phase diagram of the model can be divided into 
three regions: a superconducting phase, an insulating 
phase, and a region of phase separation. The supercon- 
ducting phase is characterized by a homogeneous density, 
leading superconducting correlations, and nonzero single- 
particle gap A = \Eo{N) - l{Eo{N + 1) + Eo{N - 1))| 
for periodic boundary conditions (PBC). Here Eo{N) is 
the ground state energy for N particles. We find this 
phase for small and moderate values of the pair hop- 
ping 1 1^1 > 1 and all fillings except n = 1/2. At ex- 
actly half-filling, an incompressible insulating phase is 
formed with exponentially decaying superconducting cor- 
relations. For large values of the pair hopping 1 
we find phase separation with the formation of particle 
clusters. In the following we concentrate on the super- 
conducting phase and check the criteria (i)-(iv). For our 
numerical analysis, we take W = —1.8 and the filling 
n = 1/3 as representative values resulting in a homoge- 
neous superconducting phase for system sizes L = 12, 24 
and L = 36 with even number of particles. 

(i) The ground state degeneracy can be studied by 
looking at the energy gap AEn{N) = En{N) - Eo{N) 
between the ground state and the n-th excited state. As 
shown in Fig. [2^, in the case of OBC, the gap between 
the ground and the first excited state A^^i^bc closes ex- 
ponentially in the system size (left panel) indicating the 
degeneracy of the ground state in the thermodynamic 
limit. This is in contrast to the case of PBC that is de- 
picted in the right panel of Fig. [2|l (blue open triangles). 
Here we find that A£^i^pbc closes linearly in the system 
size, and A£^i^pbc = A£^2,pbC7 i-c the first and second 
excited state are degenerate (blue open and closed tri- 
angles). For OBC, A£^2,OBC also closes linearly in the 
system size (red diamonds). We find that the two de- 
generate ground states in the case of OBC differ by the 
parities of the individual wires. Note that for OBC we 
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FIG. 2. a) Closing of the energy gaps with the system size 
(L = 12,24,36) for W = -1.8 and n = 1/3. For OBC, 
the gap A^i closes exponentially (left panel), in contrast to 
the polynomial closing in the case of PBC (right panel, open 
blue triangles). The energy gap AE2 closes polynomially in- 
dependent of the boundary conditions (right panel, red di- 
amonds for OBC and closed triangles for PBC). Note that 
A^2,PBC = A^i,PBC- b) Non-local fermionic correlations 
Gij on the upper wire for L = 24. c) and d) Entanglement 
spectrum for the system of the size L = 24 shows double 
degeneracy for both OBC (c) and PBC (d). 



also have A = 0. 

(ii) The intrawire single-particle correlation function 
Gij = {ouiCij) for the system of the length L = 24 is shown 
in Fig. (2)3 for the case where / = 1,2 is close to the left 
edge and j G [l^L]. We see that G/j, being exponen- 
tially small inside the wire, attains a finite value at the 
right edge showing the existence of non-local fermionic 
correlations typical for a system with MF edge states. 

(iii) Topological order (TO) manifests itself in the de- 
generacy of the entanglement spectrum (ES) j32ll34] : Let 

pA = '^i^^'^j^^ reduced density matrix of the 

system with respect to some bipartition with support on 
both wires, where p^-^^ describes a pure state of N parti- 
cles with the corresponding eigenvalues A^-^^ . In a topo- 
logical phase, the low- lying eigenvalues A^-^^ are expected 
to be doubly degenerate for each A/", for both OBC and 
PBC, as it is demonstrated in Figs. ^ (OBC) and ^ 
(PBC) for a system of the size L = 24. Moreover, the 
distributions of the low-lying eigenvalues as a function of 
A^ share the same patter in the two cases. 

(iv) The robustness of the above properties against 
static disorder is one of the key manifestations of a non- 
local topological order. We model the disorder by adding 
the term Hy = E.- V^^'^aU^ + vj^^lb^ to the Hamilto- 

nian, where V^-^^^with j = a^b are random local potentials 
equally distributed in the interval [— K, K]- We find that 
even for moderate disorder Vr = O.lt, the ground state 
remains doubly degenerate, and the system still exhibits 
the non-local correlations (Fig. pp) as well as the degen- 
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FIG. 3. Effects of imperfections on the topological order (L = 
24, n = 1/3). a) Ground state degeneracy in the presence of 
an interwire single-particle hopping H± = JZi'^y^l^i + 
b)-d) Effects of static disorder: The non-local correlations 
(b) and the degeneracy of the ES (c) indicated the topological 
state in the presence of disorder with Vr = O.lt. d) Breaking 
of the topological phase by a strong disorder with Vr — 1.5t. 



erate ES (Fig.|3^), indicating the presence of topological 
order. For strong local disorder, however, the topological 
effects disappear, as exemplified by the non-degenerate 
ES for Vr = IM in Fig. [3]i. 

Remarkably, the observed topological order and its 
consequences are also stable against a single-particle hop- 
ping H± = tya\hi^h.c. between the two wires, which 
breaks the parity of the wires and related Z2 symme- 
try. As an example, in Fig. [3^ we show the energy gap 
^i,OBC as a function of ty'. The ground state of the sys- 
tem remains quasi-degenerate (A£^ c:^ 10~^) up to values 
ty of the order of O.lt, in agreement with the prediction 
of Refs. [20, 21 . Note, however, that the dependence 
of Ai^oBC on L changes from exponential to power law 
[22 . This stability could be very important for exper- 
imental realizations of the model because the interwire 
single-particle hopping is one of the most probable im- 
perfections. 

Pair hopping with cold fermionic atoms. The key in- 
gredient of the Hamiltonian ([T]) is the interwire pair hop- 
ping with coupling W in the absence of (parity violat- 
ing) single particle tunneling. The basic idea behind an 
atomic implementation is to introduce offsets in optical 
lattices, which suppress single particle hopping by energy 
constraints, while an energy conserving pair hopping is 
allowed and mediated by interactions. 

An atomic setup illustrating these ideas is given in 
Fig. [1] while technical details and variants of the scheme 
can be found in the SL We implement the two wires 
of spinless fermions as a bipartite lattice for spinful 
fermions. Odd and even lattice sites j trap the spin t and 
^ components of the fermions with energies 62 and — ei, 
respectively, and transitions between the adjacent wells 
are induced by an external RE field or a Raman assisted 



hopping (c.f Fig. ^p). This realizes the first line of H 
in Eq. ([T]). To understand the pair hopping mechanism, 
consider the plaquette indicated in Fig. ^ by the dashed 
line. We assume an auxiliary molecular site in the cen- 
ter of the plaquette (indicated as (C) in Fig. [l]^), which 
traps both ^ and ^ atoms, and is connected to the lattice 
sites on the wire by a spin-preserving tunneling coupling 
with amplitudes tam and t^m- Pairs of atoms occupying 
the molecular site are assumed to interact via an onsite 
interaction U. In addition, we introduce spin-dependent 
lattice offsets, which are indicated by the —61,62 for the 
lattice sites on the two wires and for t and | species, re- 
spectively (Fig. [l]i-e). Such offsets can be generated as 
Zeeman shifts of the spin states, if a gradient magnetic 
field is applied perpendicular to the wire. 

Single particle hopping between the wires is suppressed 
in this setup: consider an atom, say in the upper wire 
a in lattice site 1 with spin f- Spin-preserving tunnel- 
ing is possible via the molecular site along the diago- 
nal of the plaquette (virtual processes are indicated in 
Fig. [l]i-e). It corresponds to the process tia^tm^t265 
which is suppressed by the corresponding energy offsets 
+62,0, —61. In a similar way also the tunneling of the ^ 
atom along i2a^im^iib is suppressed by energy conser- 
vation. However, for pair hopping tiai2a^1'im^1'2biib 
the overall energy will be conserved, since the two atoms 
can exchange energy via the interaction U. After adia- 
batic elimination of the intermediate sites when [/, 61/2 ^ 
tam,bm^ the resulting amplitude for the pair- hopping term 
is Vi^ tl^tl^{l/ei - 1/62)V^ (see SI). Note that the 
pair-hopping process tiaii6^i2at26 will also be allowed, 
but does not change the number of particles on the wires, 
and thus preserves atom number parity on the wires. A 
detailed description of this pair hopping dynamics includ- 
ing possible imperfections, e.g. induced by the Raman 
couplings, can be found in the SI. In Fig. [if we present a 
numerical analysis of the pair hopping dynamics, where 
(in units of tam = hm = 1) 62 = 26i = 2, = —20 (see 
SI) . Finally, note that the engineering of pair hopping has 
further applications in cold atom systems. For example, 
the pair hopping can be used as an entangling quantum 
gate, where the hopping of one particle (control) triggers 
the tunneling of a second atom (target). Further, it has 
applications in the context of quantum simulation, e.g. 
for lattice gauge theories emulation include ring-exchange 
and rishon determinant interactions j35l ESJ • 

Detection. Finally, we address the problem of detect- 
ing the emerging Majorana states in our AMO setup. 
Following the proposals of Ref. [40 , this could be done, 
e.g., by using standard AMO detection tools like time- 
of-flight imaging and spectroscopic techniques to probe 
the ground state degeneracy and the inherent non-local 
fermionic correlations. Demonstration of a non-Abelian 
statistic of the MEs, on the other hand, requires some dy- 
namical protocols resulting in the motion of MEs around 
each other. In our setup, one could think of a general- 
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ization of the ideas of Ref. [H] relying on single-site ad- 
dressing available in current experiments with ultra-cold 
atoms [m [43] . Another possibility would be an atomic 
analog of the fractional Josephson effect [16 using a prop- 
erly shaped external potential along the x-direction. 

Conclusions. In summary, we have shown that topo- 
logical states of matter with Majorana fermion edge 
states can be created in fermionic atomic ladders with- 
out any additional reservoir or p-wave interaction, but 
with only interwire pair hopping, which could provide an 
easier, complementary way for experimental realizations. 
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SI: PAIR HOPPING MODEL: LOW-ENERGY 
FIELD THEORY 

We present here some details of the low-energy field theory 
in the lattice pair-hopping model presented in the main 
text. The Hamiltonian is (for compactness, we adopt here 
a different notation Caj = a^, ctj — bj with respect to the 
main text) 

H = - mc\jc^,j+i + h.c.) + 



(4) 



and has a [/(I) ^2 symmetry conserving both the total 
number of particles and the parity of each wire. We apply 
here the standard bosonization lattice procedure in order 
to extract the low-energv field theory [HE]. Within this 
framework, we take first the continuum limit (ao being the 
lattice spacing) and introduce 

4,7 = ^ki^^ = -^'^o) + ^L,j(^ = i«o). (5) 

The bosonization identities have the following form: 



V27rao 
\/27rao 



(6) 



where R/L implies r = ±1 and ^(a/h) = ±1, and we have 
defined 



72 ' 



and the Klein operators 

{riR,-,,riL,',} = 0, 



^/2 ' 

[rir,a,rjs,h\ 



0, 



(7) 
(8) 



in order to preserve fermionic commutation relations. The 
free part (the first line) of the lattice Hamiltonian is then 
Ha + Hs, where (A = A/S) 



Ha 



VA 

2 



/ 



(9) 



with Tomonaga-Luttinger parameters Ka = 1 and sound 
velocities va — 2ta^^ sin(/ci?,Atto) with kF,s — ttu. 
The bosonized pair-hopping operator [the second line in 
Eq. Q] reads 

4W 



^ f — \j{n)cos[V87Ti9A{x) + 
'Oj^ J ao L 

+ a^^ cos[>/87r'i^A(x)] cos[aoks]{dx(ps)'^'^ , (10) 



tt{A) _ 
^ ~ (2™, 



where 7(71) is a density-dependent coefficient and we omit 
the Klein factors. Away from the strong coupling limit 
W/t ^ 1, where interactions between symmetric and 
(gapped) antisymmetric sector drive the system towards 
phase separation, the bosonized Hamiltonian H can be 
effectively split into two distinct sectors. The symmetric 
one is described by a Tomonaga-Luttinger liquid, whereas 
the antisymmetric one by a sine-Gordon Hamiltonian with 



the mass-term for the phase field i9a- In particular, at the 
Luther-Emery point i^^ = 2, it is possible to map the an- 
tisymmetric sector to the so called massive Dirac Hamil- 
tonian [4j 



(11) 



by employing standard re-fermionization techniques [T] [5] . 
Here ^s(s = L/R) are Dirac fermions with mass m ~ W, 
which are related to the bosonic operators via — 
47rrs(/Ps,A^ whcrc Vg^R/L — il ^ud the chiral 
= {cpA/V2 ± V2^a) 



fields J£s^R/L, A 



The mapping to 

Eq. (Illl) is possible due to the specific form of the interac- 



tion in the sine-Gordon Hamiltonian involving phase fluc- 
tuations 19 A- The effective Hamiltonian ( |11[ ) corresponds 
to the continuum limit of the Kitaev chain^ , and, there- 
fore, supports Majorana edge states. Since the antisym- 
metric and symmetric sectors are decoupled, the energy 
splitting between the two degenerate ground states is ex- 
ponentially small in the system size OS]. However, cer- 
tain external perturbations, such as interchain backscat- 
tering (in the presence of interactions between the wires) 
or impurities, may change the dependence of the energy 
splitting on the system size from exponential to algebraic 
(although with an extremely small prefactor in the vicin- 
ity of the Luther-Emery point, see the detailed instanton 
discussion in Ref. (3). 



SI: EFFECTIVE PAIR-HOPPING 
HAMILTONIAN - MODEL 1 

In this and the following Section, we describe two pos- 
sible implementation schemes realizing the pair hopping 
Hamiltonian discussed in the main text. 



Single Wire 

Before discussing the two couple wires and pair hopping 
between these wires we briefly describe our setup for a 
single wire. 

We implement the single wire of spinless fermions as a 
bipartite lattice of spinful fermions in an ID-optical lattice 
as indicated in Fig. [4^. Odd and even lattice sites j trap 
t and ^ components of fermionic atoms, respectively, and 
transitions between adjacent wells (and the associated spin 
flip) are induced either by an RF field or by an optical 
Raman transition. The corresponding Hamiltonian is 



j odd 



with 6^=1,2 the energies of the atomic states representing 
the spins, Q the Rabi frequency of the Raman drive and uj 
the corresponding frequency. After a transformation to a 
"rotating frame" the above Hamiltonian can be rewritten 
with the mapping aji aj {j even) and aj^ aj (j odd) 
as 

H =-t Y 4+r«J + (-^) Y 4^^' 
r—zLl;j j odd 



Here S — oj — {ei -\- €2) is the detuning , which acts as 
an offset (superlattice) for the odd lattice sites. The tun- 
nehng ampUtude in the lattice can be identified with the 
Rabi frequency of the drive, t = Q. Note that the above 
Hamiltonian corresponds to an effective model of spinless 
fermions hopping on a ID lattice. 



Two Coupled Wires 

The above realization of a ID wire of spinless fermionic 
atoms acts as a building block for two coupled wires (with 
a straightforward generalization to many coupled wires). 
The setup we have in mind is depicted in Fig. with 
the wires denoted by a and b, respectively. In addition we 
introduce auxiliary sites between the wires, which we call 
center or molecular sites in Fig. ^p. They are the basic 
ingredient in implementing pair hopping, while strongly 
suppressing unwanted processes like single particle tun- 
neling (which violates parity of atoms on the two wires). 
We extend the single wire model of the previous subsection 
assuming that the spins are placed in a spatially varying 
magnetic field with gradient perpendicular to the wires. 
This results in a spin dependent energy offset: the spins t 
and ^ on the (upper) a wire have Zeeman energies 62 and 
— ei, and the corresponding energies on the (lower) wire b 
are — ei and 62, respectively. Thus we have the following 
Hamiltonians for the two (uncoupled) wires: 

j odd j even 

r=ztl;j even 

j odd j even 

r=±l;j even 

If we choose uj — ei -\- €2, our model reduces again to a ID 
tight binding model for spinless fermions hopping on the 
wires, as discussed in the previous subsection. The Zeeman 
offsets will play a central role in the following discussion 
of pair hopping. 

As shown in Fig.|4]3, atoms can hop from the two wires to 
(auxiliary) central sites, which can be occupied by both 
t and ^ atoms. Two atoms occupying a central site will 
interact according to an onsite interaction U, effectively 
forming a "molecule" m. The corresponding Hamiltonian 
is 

H = Ho + Hi = Hoa + Hob + Hoc + Hi . 

Here Hoa and Hob are Hamiltonians for the wire as dis- 
cussed above. The Hamiltonian for the central sites c is 

Hoc^ X Ua\^al^aica^c 

c={j,j + l) 

c=(j,j + l) a = a,b 

where we adopt the notation c = + 1) for c on the 
link j^j + 1, and we have an interaction term U between 
atoms with different spin in the center of the plaquette. 
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FIG. 4. Microscopic illustration of Model I. Panel a): the sin- 
gle (a-) wire described by H is obtained by combining a species 
dependent lattice for the t and ^ species, with Raman assisted 
tunneling between the two, the small detuning thereof provid- 
ing a finite value of the S parameter. Panel b): two- wire setup 
of the Hamiltonian H\j, indicating the correspondent energy 
off-set of each site. The upper and lower wires are denoted as 
a and b, respectively. Here, in the single plaquette included 
in the dotted box, thin lines link sites connected by standard 
tunneling, and dashed curved lines link sites connected by 
Raman assisted tunneling. 



The last line is a spin-preserving hopping from the wires 
to the central site. 

Finally, Hi accounts for RF or Raman induced spin-flip 
transitions. 

Hi = - X (q' e~'''^^ a\^aic + h.c.^ 

c=(j,j + l) 
j even 

+ a|^aj+i-^,a + a|^aj+i-^,b^ e"""^* + h.c. 

j odd 

it it \ +*^* I u 

+ a|^aj+i^,a + a|^aj+i-^,b j e + h.c. . 

The first line is a spin- flip at the central site, and the last 
lines account for an RF or Raman assisted hopping from 
the wire to the central site. Note that - in contrast to the 
RF or Raman couplings along the wires - all these terms 
are off-resonant because the drive frequency cj = ei + 62 is 
detuned from from the corresponding transition frequen- 
cies ei,e2- Thus these terms average to zero, provided 
Q.' ^Q." <^ e. We will neglect them in the following dis- 
cussion. 

To illustrate the physics of pair hopping we consider now 
a single plaquette j = 1, 2 (dotted box in Fig.|4]3). We can 
write the corresponding Hamiltonian in a frame rotating 
with cj as 

Ha = e2a\^ai,a - ^ia\^a2,a - {^a\^a2ae~'''^^^ + h.c.^ 
+ e2a\^ai^b - eia\^a2b - {^a\^a2be~'^^^^ + h.c.^ 

- X tma {a\^ajc, + a|c«jc. + h.c.^ 

a.=a,h 
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For the sites on the two wires we have again used the 
notation ai^,a aia, ai^^b aib and a2i,a a2a, 
tt2t,fe ~^ CLj,b- The first line is again the Raman or RF 
hopping between sites 2^1 with energies — ei and +62, 
respectively, due to absorption (emission) of a photon 
uj = €1+62, which is tuned to compensate the energy 
difference. The second line corresponds to the central site 
with interaction and hopping. The last two lines describe 
a spin-flip on site c, and a Raman assisted hopping and 
accompanying spin flip from the wires to the central site 
c. 

Let us now analyze the various processes on the plaquette 
according to the above Hamiltonian. We will argue that 
an adiabatic elimination of the central site will result in an 
effective Hamiltonian for the two wires where single parti- 
cle tunneling is suppressed, while atoms can hop pairwise 
between the wires. 

• Suppression of single particle interwire hopping. 
Consider a single atom on the plaquette which 
occupies initially, say, site la, i.e. has spin t- 
The particle hopping to and from the central site 
preserves spin, and thus it can only tunnel along 
the diagonal 1 ^ c 1 ^- In view of the 
energy mismatch +62, 0,— ei this tunneling process 
is not energy conserving and thus will occur only 
as virtual process, which renormalizes the (single 
particle) tunneling parameters and onsite shifts of 
the wire sites. 

[Note: The above argument ignores the effect of 
the drive on the central site. We argued above 
that this coupling will be small for <^ 2e, but 
can result in a (weak) energy conserving transition 
1 'f^a 1 ^,b. We can suppress such terms by 
tilting the plaquette, so that the energies of the four 
lattice la, 2a, 16, 16 sites are +ei, — 62? +^2? with 
e'l — €2 ^ tma, This makes the spin- flip interwire 
tunneling an energy non-conserving process, i.e the 
tunneling terms of the form a^^a^b will be absent.] 

• Interwire pair hopping. While the single particle in- 
terwire hops 1 ^ 1 ^ and 2 ^, a ^ 1 6 are 
individually forbidden, the joint hopping is energet- 
ically allowed. For this pair hopping to happen the 
two particles must occupy simultaneously the central 
site to be able to exchange energy, i.e. to interact. 
An adiabatic elimination of the central site therefore 
gives a term W a\^alj^aiaa2a -\-h..c. , where in the limit 



of large U 



W = 



2 .2 .2 
u ' 



Particle assisted tunneling, where, for example, a 
particle hops 1 ^^a ^ 2 ^^h while a second particle 
2 I, a ^ 2 I, a remains on wire a, is suppressed by 
energy conservation. A process 1 t, « ^ 2 t, 6 while 
2 ^, a ^ 1 ^^a requires a spinflip on the central site 
which we argued to be small; in addition this process 
can be suppressed by tilting the plaquette. 

• Parity-preserving perturbations. Imperfections which 
preserve the parity symmetry are also present. Local 



off-sets of the form: 

jodd jeven 

+ XI ^it,ba]t,b«^t,b + XI i3i,ha]^^^aji,b (12) 

jeven jodd 

are also generated within second order perturbation 
theory, the corresponding coefficients being: 

4.2 ,2 



£2 ei 
,2 ,2 

^ ^am ^ ^bm 

i^ji,a — 5 ^j-l,b — 

€2 £2 



(13) 



Moreover, in fourth order perturbation theory, addi- 
tional diagonal interactions emerge, induced by vir- 
tual processes where two particles from the site j, j' 
belonging to the same plaquette hop into the inter- 
mediate sites, and subsequently hop back to j^j' . 
The corresponding terms read: 



cr,(T' ,^ a,a' = a,b 

with coefficients: 



f^jcr,a.'^j' a' ,a.' 



K 



1 1 \ ^ /2 .2 



€2 ei J U 



(14) 



are also present. Regarding the relevance of such ad- 
ditional terms in the effective Hamiltonian, we refer 
the reader to the discussion of the imperfections in 
the following section. 



Numerical analysis of a single plaquette 

In the following we investigate a single plaquette and 
carry out a real-time evolution under the Hamiltonian 
H\j starting from a state with two particles on the up- 
per wire, 1^) = ajaJlO). We plot, as a function of 
time, the double occupancy Pa{t) = {na,i{t)na,2{t)) and 
Pb{t) = {nb,i{t)nb,2{t)) in the upper resp. lower wire. A 
large value of Pb{t) is an indication of pair hopping. Sin- 
gle particle hopping, which is indicated via a finite value 
of S{t) = j=i 2(^«>*(^)^^.i(^)) should be suppressed. 
Taking tam = tbm = 1 as the unit of energy, we depict, 
in Fig. [5^), numerical results for Pa,b{t) and S{t), where 
£2 = 2ei = 2, [/ = -20, e' = 1/2{U - (e2 - ei)) + 6, and 
^ = 1. We find a large and finite value for the pair hop- 
ping, while the single particle hopping is of the order of 
10~^. The occupation of the central sites (Fig.jHJ)) is very 
small. 



Relation to Model 2 

As a final remark we will relate the present model to the 
one described in the following subsection, where the role of 
tunneling and Raman beams is interchanged. We can make 
a time-dependent transformation to rewrite the above pla- 
quette Hamiltonian as 
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t/tam t/tam 

FIG. 5. Real time evolution for the parameters in the main 
text, a) Time evolution of the expectation values Pa{t) = 
{na,i{t)na,2{t)) and pb{t) = (n^,! (t)n6,2(t)) that indicate the 
pair hopping. Single particle hopping indicated by S{t) — 
j=i 2(^"'^(^)^^>i(^)) order of 10~^. b) The oc- 

cupation of the intermediate sites (riac), {^hc)^ is very small 



ifn = - [p^al^aia + h.c. j - \ Q.a\^^a2h + h.c. j 
+ Ua\^a\^aica^c 

, { \ -i(e+e')t I t i(e-e')t , i \ 

- tma {a\.c^iae ' + a|^a2ae ^ ' +h.c.l 

- tmh ( aj^aifee^ ^ ' + a|^a26e ^ ^ ^ + h.c. 1 . 

The first line now looks like a tunneling Hamiltonian, while 
the last two lines correspond to time-dependent hoppings 
to the central site with frequencies e =1= e^ In this rewrit- 
ing of the model the lattice offsets have been converted 
to effective time-dependent RF or Raman couplings. In 
the following section we investigate a model, which is a 
generalization of this scheme. 



SI: EFFECTIVE PAIR-HOPPING 
HAMILTONIAN - MODEL 2 

Building block: Single wire with resonant coupling 
to intermediate states 

The basic building block of model II is represented by a 
single plaquette, illustrated in Fig. [6^, where the upper 
sites belong to the a- wire and the lower sites to the h- one, 
while the middle sites are employed as intermediate sites 
(which will then be adiabatically eliminated) to generate 
the pair tunneling term. Before discussing the emergence 
of the full pair-tunneling Hamiltonian, we show how one 
can couple one of the wires to the intermediate sites such 
that pair tunneling from the wire to the middle sites is 
generated in second order perturbation theory. This dis- 
cussion will then be naturally extended to the full pla- 
quette treatment, where the emergence of inter-wire pair 
tunneling will emerge as a combination of a pair tunneling 
process from the a-wire to the central sites, and from the 
central sites to the 6-wire. 

We start by analyzing the system composed by the lower 
part of the plaquette in Fig. [6^. Atoms in the atomic state 
I a) can occupy the 'wire' sites |aj), and the central 

site |ac); atoms in the atomic state \h) can occupy the 
central site |6c)- The microscopic Hamiltonian H for the 



subsystem can be decomposed as a sum of three terms: 
Ha^Ka + Hc + Hn,a- (15) 

The first term describes tunneling of the fermions within 
the wire: 

Ht^a = -ia{a]aj+i + h.c), (16) 

where ia is the tunneling amplitude, and a](aj) are cre- 
ation (annihilation) operators of fermions in the state \aj). 
The second term 

He = VaaUc + Vhhlhc + UaUcblbc, (17) 

describes the two intermediate states \ac) and \bc) in the 
middle of the plaquette. Here, aj(ac) and bl{bc) are cre- 
ation (annihilation) operators of fermions for the interme- 
diate states \ac) and \bc), respectively, Va and Vb are the 
corresponding potential off-sets, and U is the interparticle 
interaction. The last term 

Hn,a = Ja{al{aj + a^+i) + h.c.)+ 

- h{bl{Qijaj + Qi,^a^+i)e"^^i* + h.c.) (18) 

describes the coupling between the wires and the states in 
the middle of the plaquette: Atoms in the \a) state can 
tunnel from the wire to the intermediate site |ac), with 
a tunneling coefficient Ja, and can be transferred to the 
state \bc) via a Raman process characterized by (space- 
dependent) Rabi frequency Qij and detuning huui. We 
work here with time-dependent fields assuming a rotating 
wave approximation ^ cji). Note that the same re- 

sults can be obtained by describing Raman process as an 
auxiliary quantized single photon mode. 
We are now interested in the dynamics of such basic build- 
ing block in the regime where single particle occupation 
in the intermediate sites is suppressed as a far-off reso- 
nant state, while double occupancies of the form aj6j|0) 
are allowed (here, |0) is the fermionic vacuum). In order to 
illustrate it, we perform a quasi- degenerate perturbation 
theory (within the rotating wave approximation) in the 
limit \U\, huoi, Va,Vv ^ ta, \hQij/j^i\, with the additional 
quasi-resonant condition: 

u = -Va-Vt-nwi + 6u, \Su\-iaAhni\. (19) 

In this limit, the two states a]a]^^|0) and albl\0) are the 
only quasi-degenerate with total number of fermionic par- 
ticles equal to two. After eliminating the states with singly 
occupied intermediate sites in the second order perturba- 
tion theory, we obtain the following effective Hamiltonian: 

^eff,2 ^ -taia^jaj^i + h.c.) + Waialblajaj^i + h.c.) + 

+ 6u{alacblbc) + ^a,ja]aj + ^aj+ia'^j^^aj^i, (20) 
with a renormalized tunneling rate 

ta = ia - {Jaf/Va - ^1,, + 1 / (K - /^l), (21) 

an effective interaction in the intermediate site 5u ^ and a 
potential off-set on the a-wire: 

t _ _ l^^a^J I (cycy\ 

Va (H+fi^l)- ^ ^ 
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Moreover, a pair tunneling from the wire to the central 
sites emerges, with coefficient: 



1 



Va Vh + fvUOl 



hJai^lJ + l (23) 



Note that in order to have a finite coefficient of the pair 
tunneling, space-dependent Rabi frequencies are required. 
This is due to the fact that, for a short-range interpar- 
ticle interaction, one has to use spatially inhomogeneous 
Rabi frequencies in order to change the symmetry of the 
spatial part of the wave function for two particles from an- 
tisymmetric (for two initial particle on one of the wires) to 
symmetric (for two particles in the center of a plaquette). 
In a similar way, one describe the system of the 6-wire 
(upper part of the plaquette) coupled to the intermediate 
sites. The corresponding microscopic Hamiltonian reads: 



Hb = Ht,A,b + He + HciJ 



(24) 



where 



Ht,A,b = -ib{b]bj+i + h.c.) - A{b]bj + 6]+i6,+i), (25) 

describes tunneling in the 6-wire, and an off-set A. Here, 
b^j{bj) are creation/annihilation operators corresponding 
to the states \bj). The term 

H^,b = -h{al{Q2,jbj + Q2,j+i^j+i)e''^^* + h.c.) + 

+ J5(6j(6,+6,+i) + h.c.), (26) 

describes the tunneling between the 6-wire and the inter- 
mediate site \bc) with the amplitude Jb, and the coupling 
between the 6- wire and the intermediate site \ac) via a Ra- 
man process with detuning huj2 and Rabi frequency 0.2,j. 
The energy off-set is introduced in order to suppress single 
particle tunneling between the wires (as discussed below) . 
This off-set, however, has to match the condition 



2A = h{uji +CJ2) 



(27) 



in order to ensure resonant pair-tunneling. 
After eliminating the states with singly occupied inter- 
mediate sites in the second order perturbation theory in 
\U\,huj2,Va,Vb ^ t,hQ, we obtain the effective Hamilto- 
nian: 

i^f '2 = -tb{b]bj+i + h.c.) + Wb{alblbjbj+i + h.c.) + 

+ SuiaUcblbc) + ^b,jb]bj + 6,,+i&Li6i+i (28) 



with 



Wb 



1 



+ 



1 



Vb + A Va-nw2 + A 



hJb{^2,j-\-l — ^2,j) 



tb = k - {JbfKVb + A) - Q.2,j+i9.lj/{ya -Huj2 + A), (29) 



(14 + A) {Va + A + nw2)- ^^^^ 

Despite its simplicity, this setup already displays the fun- 
damental features of the emergence of inter- wire pair tun- 
neling. We will now show how the two wires can be res- 
onantly coupled to the intermediate sites while avoiding 
single particle inter-wire tunneling. 
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FIG. 6. Panel a): Schematic Hilbert space of a single 
plaquette, including states in the upper and lower wires 
(|aj/j+i), \bj/j^i)) and states in the center of the plaquette 
|a)c, \ b)c- Panels b,c): Schematic examples of perturbation 
theory processes generating the pair tunneling term. Below 
each perturbation step, denoted as I, II, HI and IV, the en- 
ergy of the intermediate state is indicated. In case of constant 
Qi/2,j, these two processes interfere destructively. Panel d): 
Diagrams corresponding to the illustrations in panels b and c), 
where r defined the time evolution of the perturbative steps, 
and dashed (thin) lines denote tunneling (Raman) processes. 



Coupled wires 

Our goal is to derive an effective Hamiltonian for the full 
plaquette where i) single occupancies in the intermediate 
sites are suppressed, and ii) single particle tunneling be- 
tween the wires is suppressed by energy constraints. The 
first condition is fulfilled in the regime discussed in the 
previous section; moreover, as single particle tunneling 
between the wires is always driven by a single auxiliary 
field 7 , the condition: 

|A - hwi\, \A - nw2\ > iajb, \hQi,j/j+i\, |;i^^2,,7,-+i| (31) 

guarantees that inter-wire single particle tunneling are far- 
off resonant processes, and thus suppressed within pertur- 
bation theory (see Figjs]). 
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FIG. 7. Schematic two-particle level scheme in a single plaque- 
tte (here, V^a = H > in the upper panel, and Va = Vb < 
in the lower, for the sake of clarity). The indices l,i run over 
+ 1), indicating the multiplicity of the various interme- 
diate states (15 states are involved). The Raman detunings 
huJi/2 are indicated by the green/red line, respectively, show- 
ing the two particle resonant condition in Eq. (27). 



In the perturbative regime \U\,huji,huj2, A,Va,Vv ^ 
ta,tb,\hQij/j-^i\,\hQ2,j/j+i\, after combining the results 
of the previous section, one obtains the following effective 
Hamiltonian: 



= -ta{a]aj+i + h.c.) - tb{b]bj+i + h.c.) + 

+ Wa{alblajaj+i + h.c.) + Wb{alblbjbj+i + h.c.) + 

+ 6p{alacblbc) + ^aja^aj + ^a,j+ia]^-^aj+i + 

+ ^b,jb]bj+^b,j+ib]^,bj^i. (32) 



Here, the first line describes renormalized intrawire tun- 
neling; the second one contains terms coupling the a and 
6- wire to the intermediate sites, respectively. Finally, in 
analogy with the previous discussion, the last two lines 
describe a residual interaction in the intermediate sites 
due to the quasi-resonant condition in Eq. (19), and lo- 



cal potential offsets. Beyond constituting the basis for the 
realization of the pair-tunneling Hamiltonian, Eq. ( 32 ) dis- 
plays interesting physics by itself, as its effective descrip- 
tion in terms of two Luttinger liquids coupled with a ID 
superconductor (as an effective description of the interme- 
diate wire composed by the center sites of each plaquette) 
has been also investigated in the context of emergent Ma- 
jorana edge states 0. 
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FIG. 8. Single particle level scheme, illustrating how the off- 
resonant condition between \a)j and \b)j can be engineered 
by fulfilling Eq. (IsTl). 



Inter-wire pair tunneling 

Adiabatic elimination of the intermediate sites 

The pair tunneling Hamiltonian discussed in the text is 
then obtained by adiabatically eliminating the central sites 
from Eq. (32 ) in the regime 6U ^ t, HQ, where the quasi- 
resonant condition in Eq. ( |19| ) is not met anymore. The 
effective Hamiltonian will then read: 

^eff,4 ^ -ta(a]aj+i + h.c.) - th{b]bj+i + h.c.) + 
+ W{b]b]^^ajaj+i+h.c.)+ 

+ ^,,jb]bj + + H'^^, (33) 

where the second line describe the pair tunneling term 
between the wires (which, in perturbation theory, emerges 
as a sum of terms similar to the ones illustrated in Fig.[6]), 
with coefficient: 



W 



+ 



Vb + huji 
1 



^h'jaJb 
6U 



+ 



1 



A + Va-nw2 A + Vb 
X (^i,j+i — Qi,j){Q2,j — ^2,^+1)5 



(34) 



while denotes additional imperfections (such as 

intra-wire and inter-wire interactions) emerging from ad- 
ditional fourth order processes. While all of them preserve 
parity symmetry, some of them may change the energy 
splitting between the two degenerate states in the topo- 
logical region discussed in the text from exponential to 
power law as a function of the system size. 



Additional terms 

Let us now discuss the imperfections described by H^^^ . 
Intra-wire interactions. These interactions are gener- 
ated via virtual processes, where, e.g., two particles on 
nearest-neighbor sites in the a- wire are transferred to the 
intermediate sites \a)c and \b)c via tunneling and Raman 
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assisted process, respectively, and then move back to the 
same two sites. The effective contribution in the Hamilto- 
nian reads: 



GaTlajnaJ + l + GbUbjUbJ-^-l 



(35) 



and the corresponding pre-factors in fourth-order pertur- 
bation theory are: 



Ga 

and 
Gt 



1 



Va 



+ 



1 



X 



(36) 



A + 14 



+ 



{hJt 



X (|Q2,,f + |n2,, + l|' 



huj2 + K + A 



6U 



(37) 



■ 1^2, j)- 



The effect of these terms is to renormaUzed the single wire 
Luttinger parameter, producing only some quantitative 
shifts in the model Hamiltonian phase diagram. In case 
Gb,Ga < 0, the topological phase is expected to emerge 
at smaller values of VK, as the value of Ka becomes larger 
when an additional intra-wire attraction is introduced 
Inter-wire interactions. Inter- wire terms emerge as well 
in fourth order perturbation theory. They are induced by 
virtual processes, where a particle hops from \a)i to \a)c 
and another one hops from \b)i to \b)c (or similar pro- 
cesses induced by the Raman couplings), and then both of 
them hop back to the original wires in \a)k and \b)r. The 
contribution in the effective Hamiltonian then reads: 



{Ki^i^k,ralaiblbr -\-h.c.). 

i,l,k,r—j,j-\-l 



(38) 



Here, the indices i,l,k,r belong to the same plaquette, 
that is, i, l,k,r ^ {j, j The corresponding coefficients 
are: 



l,k,r 



1 



+ 



1 



Va Vt + A 
1 



t2 t2 



Vb + Va + A + U 
2 



Va -huJ2 + A Vb + huJl 
^^^l,i^2,r^l 2^2 k 

Vb + Va + A + u + nwi -huJ2 



X 



+ 



(39) 



While not breaking the parity symmetry, these terms may 
lift the splitting of the ground and first excited state from 
exponential to algebraic as a function of the system size [3^ ; 
as such, it'd be preferable to strongly reduce their effects 
in order to guarantee observability of the topological fea- 
tures even in small size systems. In general, the ration 
between W and Ki^i^k,r can be minimized by tuning the 
Hamiltonian parameter in such a way that the sum of all 
processes involved in Ki^i^k,r is annihilated via quantum 
interference; this requires matching conditions involving 
both interactions, off-sets and detunings. A simpler, al- 
ternative way is to tune the potential offsets such that 
Va = —A — Vb, so that the first contribution is killed via 
interference, and then take J ^ 4Q; this way, the ratio 
W/Ki,i,k,r will be of order 1/(4 * (16)) = 1/64, thus 
negligible with respect to other possible imperfections in 
the system. 
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FIG. 9. Real time evolution for the parameters in the main 
text and 6U = OkHz, SV = 0.02kHz. a) Time evolution of 
the expectation values Pa{t) = {na,i{t)na,2{t)) and pb{t) = 
{nb,i(t)nb,2(t)) that indicate the pair hopping. Single particle 
hopping indicated by S{t) = J2i j=i 2(^«'^(^)^^>i (^)) 
order of 10"^ b) Since the detuning ^ = 0, the occupation of 
the intermediate sites (riac), (^bc)? is relatively large large. 



Other contributions. Within the effective Hamiltonian, 
the additional potential off-sets should have the same ef- 
fects as the intra-wire terms, and it can be shown that 
they strength can be tuned by slightly shifting the ini- 
tial potential off-set A from exact degeneracy according 
to Eq. (27). Moreover, it can be shown (by employing a 
fioquet formalism) that imperfect addressing of the indi- 
vidual sites by the Raman beams does not qualitatively 
modify the effective Hamiltonian, as the quasi-degenerate 
subspace within perturbation theory is not affected. 



Numerical results 

In the following we present quantitative results on model 
I defined on a single plaquette. Starting from a state with 
two particles on the upper wire, |^) = ajaJlO), we carry 
out a real time evolution under the Hamiltonian H. Based 
on our numerics, we answer the following two questions: 
(1) Can we suppress single particle hopping between the 
two wires? (2) Do we get pair hopping? Further, we check 
the validity of our perturbation theory assumptions. First, 
we investigate the effect of giving up the space dependence 
of the Rabi frequencies and take = ^j+i- We expect 
that the pair hopping disappears. Second, we check the 
resonance condition for pair hopping, A = ^{^i +^2), 
by comparing to a protocol where A = huoi, i.e. sin- 
gle particle hopping is resonant. For the numerical inves- 
tigation, we use the experimentally realistic parameters 
ta = tb = lOHz, Ja = Jb = = IkHz, Vb -6kHz, 
huoi = 20kHz, h(jj2 = 12kHz. The resonance conditions 
imply that 

Va = -A-Vb + 6V, (40) 

U = -hwi -Va-Vb + SU, (41) 
where we have defined SU and AV to account for small 
deviations from the resonance conditions, as the are ex- 
pected in a realistic experimental setup. 



Suppression of single particle hopping 

We carry out a numerical analysis that proves the sup- 
pression of single particle hopping in model I using the 
parameters defined above. As a figure of merit, we take 
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FIG. 10. Real time evolution for the parameters in text 
and 5U = 0.5kHz, SV = 0.02kHz. a) Time evolution of 
the expectation values Pa{t) = (na,i(t)na,2{t)) and Pb{t) = 
{nb,i{t)nb,2{t)) that indicate the pair hopping. Single parti- 
cle hopping indicated by S{t) = 2(^a,*(^)^b,j(^)) 
the order of 10~^. b) Occupation of the intermediate sites, 
(riac), (^6c) a function of time. The detuning S leads 
to relatively small occupation of the intermediate site. 
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FIG. 11. Real time evolution for the parameters in the main 
text and single particle resonance, A = huji. a) Time evo- 
lution of the expectation values Pa{t) = {na,i{t)na,2{t)) and 
Pb{t) = {nb,i{t)nb,2{t)) that indicate a vanishing pair hopping, 
b) Single particle hopping indicated by S{t) is large. 



FINAL REMARKS 



the expectation value S{t) = j=i 2{'^a,i{t)nb,j{t)) that 
indicates if we find, for some time t, one particle on the 
upper and the other particle on the lower wire. We find, 
that for 6U and 6V on the order of a few hundred Hertz 
the expectation of S{t) is of the order of 10~^. 



Creation of pair hopping 

Next, we present a numerical proof that we can engi- 
neer pair hopping in our setup. As a figure of merit, 
we take the expectations Pa{t) — {na,i{t)na,2{t)) and 
Pb{t) = {nb,i{t)nb,2{t)) indicating that the two particles 
are both on the same wire. If pair hopping is possible, then 
Pa{t) should decrease with time, while Pb{t) increases. In 
Fig. [9] we present results for SU = 0, SV = 0.02kHz. We 
see that the occupation of the intermediate site is rela- 
tively large. This problem can be overcome by introduc- 
ing a detuning SU / 0. As shown in Fig. [To] already a 
small detuning of SU = 0.3kHz, SV = 0.02kHz signif- 
icantly reduces the occupation of the intermediate sites. 
Further, from the width of Pa{t) and Pb{t) we can deduce 
that W - 1/lOOkHz = lOHz is of the order of the hop- 
ping ta. 



Check of the perturbation theory 

Now, we investigate the validity of our perturbative ap- 
proach. First, we give up the space dependence of the Rabi 
frequencies and take Qj = ^j+i- We expect that the pair 
hopping disappears, and that there is no hopping between 
the two wires. This assumption is confirmed by our nu- 
merical analysis: The expectation Pa(t) — 1 for all times, 
while pb{t) = 0. Further, S{t) is of the order of 10"^ 
Finally, we consider the case where we tune the parameters 
in our model such that they are resonant to single particle 
hopping, i.e. we take A = huJi. We find that Pb{t) = 
for all times in agreement with a vanishing pair hopping 
(see Fig.pT^i)). The results for the single particle hopping, 
indicated by a non- vanishing values of S{t) are depicted 
in Fig.[TT]3). 



Let us conclude with some final remarks on the realiza- 
tion of the pair hopping Hamiltonian. As discussed in the 
previous sections, pair tunnelings of order W/t :^ 1 can be 
achieved when considering tunneling rates tjh 10 — 20 
Hz (as currently employed in ultracold fermion experi- 
ments [8^) for Rabi frequencies of order Q ^ 27r x 6 kHz. 
Employing Raman beams with larger Rabi frequencies 
would lead to even better energy scales (although suffi- 
ciently strong interactions would be needed as well). An 
alternative route toward the observation of Majorana edge 
states would be to consider the plaquette Hamiltonians 
without tracing out intermediate states with pair of par- 
ticles in the intermediate sites, that is, obtaining Hamil- 
tonian of the form in Eq. ( |32| ), where the corresponding 
couplings Wa, Wb may be much larger that W , as they are 
derived in second order perturbation theory. The numeri- 
cal results confirm this enhancement (in the specific case, 
by a factor of ^ 10). The same considerations apply to 
model II. The corresponding many-body scenario would 
then be very similar to the one described in Ref. 3, where 
the a and 6- wire are couple to an additional ID s-wave su- 
perfluid with a large spin gap (described by the dynamics 
in the intermediate sites). 

Minimal instances of the present scheme can be validated 
on a single plaquette, by measuring the relative parity of 
a certain state under the evolution of the effective Hamil- 
tonian by using band-mapping or in-situ imaging tech- 
niques 0. Moreover, a generalization to multi-leg ladders 
and 2D setups is indeed conceivable under the same as- 
sumptions described in the previous sections. 
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